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!>• ! Abstract 



The dilepton production at backward rapidities in pAu and pp collisions at RHIC and LHC ener- 
gies is investigated in the dipole approach. The results are shown through the nuclear modification 
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. ratio R p A considering transverse momentum and rapidity spectra. The dilepton modification ratio 

t> ■ 

' presents interesting behavior at the backward rapidities when compared with the already known 

forward ones, since it is related with the large x kinematical region that is being probed. The 



rapidity dependence of the nuclear modification ratio in the dilepton production is strongly depen- 



Qh! dent on the Bjorken x behavior of the nuclear structure function ratio Rp 2 = F^/F^- The RpA 

4^ ■ transverse momentum dependence at backward rapidities is modified due to the large x nuclear 



effects: at RHIC energies, for instance, the ratio R p a is reduced as pr increases, presenting an 
opposite behavior when compared with the forward one. It implies that the dilepton production 
at backward rapidities should carry information of the nuclear effects at large Bjorken x, as well 
as that it is useful to investigate the px dependence of the observables in this kinematical regime. 
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I. INTRODUCTION 

The hadron transverse momentum spectrum measured by the RHIC experiments [lj] is one 
of the most accessible distributions carrying information about the high density nuclear envi- 
ronment. At forward rapidities, the transverse momentum distribution for hadrons presents 
results compatible with the Color Glass Condensate (CGC) description of the saturated 
regime at high energies [2|. These results are investigated through a nuclear modification 
ratio, which relates dAu and pp collisions. It is well known that the transverse momen- 
tum dependence of this ratio shows a pronounced Cronin peak at central rapidity, which 
is suppressed at forward ones. For central rapidities the Cronin peak is due to the multi- 
ple scatterings of the projectile constituents with the dense target. At large rapidities the 
Cronin peak suppression is understood as a signal of the saturation phenomena intrinsic of 
the Color Glass Condensate (CGC). However, at backward rapidities, the recently measured 
hadron nuclear modification ratio shows a pronounced peak |3j, which should be due to final 
states interactions. 

In this work the dilepton production at high energies in the backward rapidity region 
is investigated in the proton-nucleus and proton-proton collisions. The analysis at forward 
rapidities was already done 4j and the results show that the dilepton is also a suitable 
observable to investigate the Color Glass Condensate. Here, the interest to explore the 
backward region in the dilepton sector relies on the study of the nuclear effects at large and 
small x and on the comparison with the Cronin effect in this region. In this kinematical 
region, the nucleus interacts by means of large Bjorken x partons, and no saturation effects 
are expected consequently. In this regime, the proton interacts through the small x partons. 

It has been shown (4] that the dilepton production at forward rapidities, analyzed in the 
context of the Color Glass Condensate in pAu collisions, presents the same features of the 
Cronin effect at forward rapidities, implying that this should be considered as an initial state 
effect. However, at backward rapidities, the Cronin peak present in the hadrons RHIC data 
{3] is still not understood, and could be related to final state effects. The dileptons analyzed 
in this kinematical region, since they do not interact strongly with the final environment, may 
not present the same behavior of the hadron spectra. Although it could imply the dilepton 
production to be symmetric when forward-backward regions are compared, in contrast to 
the different behavior of the hadron spectra in both regions, we verify that the dilepton 
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production at backward rapidities is strongly modified due to large x nuclear effects. 

This work is organized as follows. In the next section the dipole approach is presented 
for the dilepton production at backward rapidities. In the Sec. 3 the nuclear partonic 
parametrizations employed in this work as well as a comparison between then are presented. 
In the Sec. 4 the results are discussed. Our conclusions are left to the last section of this 
work. 

II. DILEPTON PRODUCTION AT BACKWARD RAPIDITIES 

To evaluate the dilepton production at backward rapidities, an adequate treatment of 
the approach employed for the forward ones needs to be considered. The investigation at 
forward rapidities in pA collisions requires to describe the nucleus as a high density saturated 
system (CGC), and the dileptons are produced by the decay of a virtual photon emitted 
by the quark from the proton, which interacts with the nuclear dense medium (basically a 
dipole approach in the momentum space with a specific configuration of the nucleus). In 
this case the dilepton cross section can be factorized considering the wave function of the 
fluctuation q / ~f* and the interaction cross section of a quark with the nucleus, which takes into 
account the high density effects. However, at backward rapidities the nucleus interacts by 
means of large Bjorken x partons, being no more described by a high density QCD system. 
The proton interacts through very small x partons implying that the dipole picture could 
be applied in this kinematical regime. 



Nucleus 




FIG. 1: In the dipole approach, dilepton production looks like a bremsstrahlung. A quark or an 
anti-quark from the projectile hadron (nucleus) scatters off the target color field and radiates a 
photon (7*) with mass M (before or after the quark scatters), which subsequently decays into the 
lepton pair. 
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In the dipole approach at backward rapidities, the quark of the incident hadron (nucleus) 
fluctuates into a state containing a quark and a massive photon. The interaction with the 
target (proton) provides that the virtual photon is freed, decaying into a lepton pair (sf] (one 
of the diagrams is shown in the Fig. [T]) . 

We are interested in the study of the dilepton production in a region where the interaction 
time of the projectile quark with the proton target is much shorter than the q^* fluctuation 
time, meaning large coherence length. In order to apply the dipole picture at backward 
rapidities, we need to exchange the nucleus with the proton when comparing with the dipole 
picture at forward rapidities. In doing such procedure, the q^* fluctuation time will be 
larger than the interaction time, since the coherence length l c has the dependence l c ocl/x\ 
in the case studied here. Indeed, the coherence length is important in the study of nuclear 
targets, where it is a fundamental quantity controlling nuclear effects jfj]. For the nuclear 
targets large coherence length needs to be reached if shadowing effect are under investigation, 
which implies that the approach should be applied only to smaller Bjorken x Q]. Another 
interesting point that arises here is that the pt broadening due to the nuclear multiple 
scattering verified at forward rapidities jf| should not be present at backward rapidities, 
since the nuclear projectile will interact with a proton target and will not probe a large and 
very dense system. 

The advantage of this formalism is that the dilepton cross section can be written in terms 
of the same color dipole cross section as small-x Deep Inelastic Scattering (DIS). Although 
diagrammatically no dipole is present in bremsstrahlung, the dipole cross section arises 
from the interference of the two bremsstrahlung diagrams, as shown in Ref. 8] in a detailed 
derivation. 

The cross section for the radiation of a virtual photon from a quark (with momentum 
fraction X2) of the nucleus scattering off a proton at high density at backward region can be 
written in a factorized form as [j], [h| , 

da DYback a 2 eri 
dM 2 dyd 2 p T ~Q^M 2 

where px is the dilepton transverse momentum, M is the dilepton mass, y the rapidity, p 



POO 

/ dpW(x2,p,Pr)<7dip(xi,p), (!) 
Jo 



is the dipole transverse separation. The variables x\ and £2 are defined in the usual way 



X(i 

\2 



— ~^ Pt e ±y ', with s being the squared center of mass energy. As can be seen from 
the definition of x\ and 22, backward rapidities will provide large X2 and small X\. We 
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are investigating dilepton production in pp and pA collisions. In a symmetric collision, the 
lab frame is equivalent to the center of mass frame considering colliders. In the case of 
asymmetric collisions, e.g. pA, the center of mass moves longitudinally in the lab frame. It 
provides that a largest interval in Bjorken x could be explored in the asymmetric collisions 



ll| if the lab frame is under consideration. In this work the ratio between pA and pp 
collisions will be investigated considering the rapidity in the lab frame. 
Based on the dipole approach, the function W(x 2 ,p, Pt) is given by [9|, 



j X2 Ot ol y 
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In the above equation a is the momentum fraction of the quark carried by the virtual 



photon, rf = (1 
Tj are given by, 
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Here, a nuclear structure function F 2 A {^- ) M 2 ) with a x 2 dependence is considered in order 
to take into account the nuclear interaction. For pp collisions the nuclear structure function 
F 2 A (^-, M 2 ) needs to be replaced by the proton structure function F^i^-, M 2 ). The dipole 
cross section is evaluated with the argument x\ and the same dipole cross section extracted 
from the deep inelastic scattering HERA data should be applied. 

The dipole cross section in the phenomenological model proposed by Golec-Biernat and 
Wusthoff 12] (GBW), which has produced a good description of HERA data in both inclu- 
sive and diffractive processes, was employed here. It is constructed interpolating the color 
transparency behavior a dip ~ p 2 for small dipole sizes and a flat (saturated) behavior for 
large dipole sizes a dip ~ o"o (confinement). The expression has the eikonal-like form, 

p 2 QI 



1 — cxp 



4(xi/x ) ; 



(3) 



where Qq = 1 GeV 2 and the three fitted parameters are <7o = 23.03 mb, xo = 3.04 x 10' 
and A = 0.288. 



5 



012345012345 

P T (GeV) p T (GeV) 



FIG. 2: Dilepton transverse momentum distribution for pp and pA (Au - considering EKS 
parametrization) collision at RHIC energies for dilepton mass M = 6 GeV. 

The pp and pAu transverse momentum distributions at RHIC energies for backward ra- 
pidities can now be evaluated, considering a specific parametrization to the nuclear structure 
function F^(x, M 2 ) in the pA case. Two distinct parametrizations for the nuclear structure 
function will be used, which will be properly discussed in the next section. In this section 



EKS nuclear parametrization 13| was employed together with the GRV98 parton distribu- 
tion function 14J in order to obtain the px distribution for the pA collision and the F%(x, M 2 ) 
is obtained from the GRV98 parametrization for the pp collision. In the Fig. [2] the pp and 
pA transverse momentum distributions for RHIC energies are shown and significative visible 
differences between both distributions are not found. For this reason, the ratio between pA 
and pp cross section should be useful to investigate modifications in the nuclear cross section 
in comparison with the pp distribution. The ratio is defined in the following way 

= dajpA) I da(pp) 
pA dp 2 T dydM/ dp 2 T dydM' 1 1 

The ratio R p a is normalized by 1/A, since the nuclear parton distributions are normalized by 
this factor. The difference between pA and pp calculations are due to the nuclear structure 
function F^, then any effect in the pt or rapidity spectra should be related to the nuclear 
effects contained in this function. 

The next-to-leading order calculation (NLO) in the collinear factorization could be eval- 
uated and compared with the results in the dipole approach. However, the evaluation of 
the transverse momentum distribution in the collinear factorization is not a straightforward 
caption and does not describe the no—ion of the Pt DteU-Yan d.leptons fl Q. 



18) for 



A resummation of large logarithms in px/M 15|, [l6| or the introduction of an intrinsic 
transverse momentum [l7| should be necessary to avoid the divergence of the differential 
cross section at px « M. New approaches have been suggested in order to describe the pt 
spectra of the dilepton (for a good discussion see Ref. 18|), however, there is a strong depen- 
dence on phenomenological parameters. Considering, for instance, the intrinsic transverse 
momentum, the px distribution is dependent on a phenomenological parameter (see Re 
for a comparison between the NLO calculations and the dipole approach, and Ref 
a new approach). Of course, the ratio R p a should be finite since the pt divergence will be 
canceled, however, a more reliable result needs to consider finite spectra. This indicates that 
the dipole approach should be more adequated to investigate the low pt dileptons at high 
energies, since there are no free parameters and the cross section is finite at low px- 

For sake of completeness, it should be stressed that some limitations of the dipole ap- 
proach at backward rapidities are in order. Here, the nucleus is considered by means of 
an integrated gluon distribution, however, at LHC energies and more central rapidities, the 
Bjorken 22 (nucleus) reaches values around 0.002, where the consideration of an integrated 
parton distribution could be questionable. At this high energy limit, the effects of finite 
transverse momentum of the incoming partons become important. A more robust treatment 
would need to include the transverse momentum of the partons from the nucleus in the ini- 
tial state of the interaction. This should be done considering the kx factorization approach 
fl9|. where the off-shell partonic cross sections are convoluted with kx unintegrated parton 
densities f a (x,k T , (i 2 )- Considering dilepton production, the kx factorization is investigated 
in the Ref. 20J and compared with a phenomenological intrinsic kx approach (kx factoriza- 
tion with on-shell partons), in order to describe the transverse momentum distribution of 
the Drell-Yan dilepton production. In spite of a reasonable data description, the kx factor- 
ization overestimates the data and the intrinsic kx approach depends on phenomenological 
parameters (two parameters). 

For the reasons presented above, we have focoused our analysis at backward rapidities, 
and not at more central ones. The use of the dipole approach together with the collinear 
factorization is justified as the goal of this work is to investigate the nuclear effects in the 
nuclear modification ratio R p a for dileptons, using nuclear parametrizations to describe the 
F£ ■ In the next section the parametrizations for the nuclear structure function will be 
discussed and compared. 
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III. NUCLEAR PARTON DISTRIBUTION FUNCTION 



In this work the parametrizations for the nuclear parton distribution function (nPDF) 
proposed by Eskola, Kolhinen and Salgado (EKS parametrization) [l^] and by D. de Florian 



and R. Sassot (nDS parametrization) [2l| are considered in order to describe the nuclear 
effects in the F^. In this section the main features of each parametrization will be discussed. 

Both nuclear parton distribution functions provide a global fitting to fixed target ex- 
perimental results, and consider the DGLAP equations for the Q 2 evolution. The initial 
conditions are adjusted to describe the DIS in lepton-nucleus collisions and the dilepton 
production in proton-nucleus collisions. Both nPDFs are used to obtain an appropriated 
parametrization for the ratio RpAx, Ql) = F^/AF^, with 

F 2 ,p = Y,<l x f^Q 2 o) + x fi ,p ^Q 2 o^ ( 5 ) 

in which f^(x,Ql) are the free proton parton distributions and fy- (x, Ql) are the nuclear 
distributions of parton flavor q. The parametrization for the parton distributions are con- 
strained by charge, baryon number, and momentum conservation. 

However, the method used by EKS and nDS differs with respect to their approaches to 
parametrize the nPDFs. In the EKS approach, the nuclear parton distributions are the 
free proton ones times a correction factor: f£ (x, Ql) = Rf(x,Ql)f^(x,Ql), and all nuclear 
effects are enclosed in Rf(x, Ql). One consequence of this definition is that the nPDFs are 
null for x > 1, although they should be non-zero for x < A. On the other hand, nDS uses a 
convolution to relate nPDFs to the free proton PDFs: 

ft^ QD = / A -W q (y, A)f> (-, Ql) , (6) 

Jx y \y j 

in which Wi(y, A) contains the information about the nuclear effects. For instance, if nuclear 
effects are disconsidered, W q (y,A) = AS(1 — y). 

The nuclear effects are verified by comparison of the nuclear structure function with 
the proton (or deuterium) structure function. In the Fig. [3] the prediction for the ratio 
F-fAF^ is presented for both parametrizations. The figure can be divided in four regions 



of Bjorken-x 22j. In the Fermi motion region x > 0.8, Rp 2 is greater than 1 and increases 
with x. The EMC region 0.3 < x < 0.8 is characterized by Rp 2 < 1. The antishadowing 
region 0.1 < x < 0.3, and the shadowing region x < 0.1 are defined by Rp 2 > 1, and 
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FIG. 3: Comparison between EKS and nDS parametrizations for the ratio F^/F^ (A=Au) for a 
dilepton mass (scale) M = 6 GeV as a function of the Bjorken x of the nucleus. 

Rp 2 < 1, respectively. The Fig. [3] shows that the parametrizations differ most in both 
ECM and shadowing regions, with EKS parametrization getting a lower ratio than nDS. 
The parametrizations are considered at the leading-order, since the approach employed here 
considers LO diagrams, without higher orders calculations. Indeed, the dipole cross section 
should carry some information about higher orders, but these properties are included in 
the phenomenological parametrization of the dipole cross section, in our case, the GBW 



parametrization 



In the next section, the results for dilepton production will be explored, considering both 
parametrizations for the nuclear parton distribution presented in this section. 



IV. RESULTS AND DISCUSSIONS 



The results for transverse momentum and rapidity distributions for the dilepton produc- 
tion (with mass M = 6 GeV) at RHIC {y/s = 200 GeV) and LHC (y/s = 8.8 TeV) energies 
at backward rapidities will be discussed. 

The ratio R p a was evaluated at RHIC and LHC energies, and in the Figs. H] and [5] the 
results are shown in 3D plots for rapidity and spectra, considering the EKS and nDS 
parametrizations. The behavior of the ratio R p a reflects the Xi dependence of the ratio 
F^/F^, presented in the Fig. [3j 

At RHIC energies a weak dependence of the ratio on the transverse momentum is verified, 
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RHIC nDS 



FIG. 4: Ratio R p a for RHIC energies considering the EKS and nDS parametrizations. 

and in general the ratio reaches smaller values at large p?, being this more evident with the 
EKS parametrization. Concerning the rapidity spectra, the EKS parametrization predicts 
a strong suppression of the ratio at very backward rapidities and large p?, in comparison 
with the nDS parametrization, which predicts an almost flat behavior. 

To explain such results it is necessary to know the range of Xi values in the rapidity and 
Pt spectra at RHIC. For RHIC energies considering rapidities from -1 to -2.6, and pr from 
1 to 7, the X2 range is between 0.08 and 0.5, respectively, meaning that for more backward 
rapidities, partons with larger X2 are being probed. The nuclear effects that appear in the 
F2 at this X2 range are mainly due to EMC effect (reduction of the ratio Rp 2 as x^ increases, 
see Fig. [3]), which provides the reduction of the ratio R p a at lower rapidities in the Fig. 0J 
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Concerning the px spectra, %2 increases with px, and as the region probed here is related 
to the EMC effect, the result is a reduction of the ratio R p a as pt increases. The large 
suppression of the ratio R p a of the EKS prediction in comparison with the nDS in the Fig. 
H] is a consequence of the large difference in the Rp 2 predictions of the parametrizations in 
the EMC effect region. 




FIG. 5: Ratio R p a for LHC energies considering the EKS and nDS parametrizations. 

At LHC energies (Fig. [5]) the rapidity spectra present a peak at intermediated rapidi- 
ties, which is more pronounced in the EKS parametrization in comparison with the nDS 
parametrization. The R p a transverse momentum dependence presents two distinct behav- 
iors: for very backward rapidities the ratio decreases as px increases and for more central 
rapidities the ratio R p a increases with px- 
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For LHC energies, rapidities from -1 to -6, and px from 1 to 7, the x 2 range is between 
0.002 and 0.3, respectively. Here we verify that not only large x of the nucleus is been probed, 
but small x too. The x 2 range probed at LHC provides that shadowing and antishadowing 
effects are present. The peak at intermediate rapidities is related to the antishadowing effect 
and the suppression at more central rapidities is related to the shadowing effect. The pt 
spectra is more involved since the ratio R p a presents different behaviors. For more backward 
rapidities the ratio is reduced for large pt {x 2 in the antishadowing region, near to EMC 
region). For more central rapidities, the ratio increases for large p?, since the x 2 is in the 
shadowing region. 

Comparing the predictions from the EKS and nDS parametrizations, in the Fig. [3] the 
EKS parametrization predicts more pronounced antishadowing, which explains the results 
found here for the ratio R p a at LHC energies (Fig. [5]). Our results presented for LHC 
energies show that dilepton carries information about the large and small x regime of the 
nuclear environment. The distinct nuclear effects are verified in rapidity and pt spectra for 
the dilepton production at backward rapidities, for RHIC and LHC energies, following the 
x region involved. 
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FIG. 6: Comparing the ratio R p a for dileptons at backward and forward rapidities for RHIC ener- 
gies, considering the EKS and nDS nuclear parametrizations at backward, and the CGC predictions 
at forward 



In order to compare the results presented in this work, with the previous results from 



the forward region 



4j, in the Fig. [6] the ratio R p a for RHIC energies is shown for positive 
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and negative rapidities. The ratio presents different behaviors concerning the transverse 
momentum dependence. While at forward rapidities the saturation phenomena implies that 
the ratio R p a increases at large p?, at backward rapidities the ratio R p a decreases as large 
Pt is reached. This distinct behavior is associated with the large x nuclear effects present 
at backward rapidities. 

V. CONCLUSIONS 

The nuclear modification ratio R p a for the dilepton production was investigated for p? and 
backward rapidity spectra in the dipole picture for RHIC and LHC energies. We have verified 
a strong dependence of the nuclear modification ratio with the nuclear effects. The results 
presented in this work are explained by the dependence of the nuclear structure function 
ratio Rp 2 on the Bjorken x. The hadron production was investigated at backward rapidities 
by PHENIX collaboration at RHIC [s| presenting the data on the nuclear modification 
ratio. The results present an enhancement in the nuclear modification ratio for 1.5 < 
Pt < 4.0, which still demands some caution, once there are some data uncertainties and a 
discrepancy between analysis methods of the data. The results presented in this work for 
dileptons indicate that the enhacement in the hadron spectra at backward rapidities could 
be mainly due to final states, since, the intial state effects investigated here do not provide 
such behavior. Moreover, the dilepton at backward rapidities is a suitable observable to 
understand and quantify the nuclear effects at large and small Bjorken x. Additionally, the 
transverse momentum dependence of the ratio R p a is strongly modified at RHIC energies, 
if forward and backward rapidities are compared, due to distinct x regions being probed. 
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